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Abstract —In this paper, we propose a new probability distri¬ 
bution function which accurately describes turbulence-induced 
fading under a wide range of turbulence conditions. The proposed 
model, termed Double Generalized Gamma (Double GG), is 
based on a doubly stochastic theory of scintillation and developed 
via the product of two Generalized Gamma (GG) distributions. 
The proposed Double GG distribution generalizes many existing 
turbulence channel models and provides an excellent fit to the 
published plane and spherical waves simulation data. Using this 
new statistical channel model, we derive closed form expressions 
for the outage probability and the average bit error as well as 
corresponding asymptotic expressions of free-space optical com¬ 
munication systems over turbulence channels. We demonstrate 
that our derived expressions cover many existing results in the 
literature earlier reported for Gamma-Gamma, Double-Weibull 
and K channels as special cases. 

Index Terms —Free-space optical systems, fading channels. 
Double GG distribution, propagation, irradiance, optical wireless, 
spatial diversity. 


1. Introduction 

F REE-SPACE optical (ESO) communication enables wire¬ 
less connectivity through atmosphere using laser trans¬ 
mitters at infrared bands. These systems provide high data 
rates comparable to fiber optics while they offer much more 
fiexibility in (re)deployment. Since they operate in unregulated 
spectrum, no licensing fee is required making them also a 
cost-effective solution m-ni. With their unique features and 
advantages, ESO systems have attracted attention initially as 
a last mile solution and can be used in a wide array of 
applications including cellular backhaul, inter-building con¬ 
nections in enterprise/campus environments, video surveil¬ 
lance/monitoring, fiber back-up, redundant link in disaster 
recovery and relief efforts among others. 

A major performance limiting factor in ESO systems is 
atmospheric turbulence-induced fading (also called as scintil¬ 
lation) jH . Inhomogenities in the temperature and the pressure 
of the atmosphere result in variations of the refractive index 
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and cause atmospheric turbulence. This manifests itself as ran¬ 
dom fiuctuations in the received signal and severely degrades 
the ESO system performance particularly over long ranges. 

In the literature, several statistical models have been pro¬ 
posed in an effort to model this random phenomenon. His¬ 
torically, log-normal distribution has been the most widely 
used model for the probability density function (pdf) of the 
random irradiance over atmospheric channels O-Q. This 
pdf model is however only applicable to weak turbulence 
conditions. As the strength of turbulence increases, lognormal 
statistics exhibit large deviations compared to experimental 
data. Moreover, lognormal pdf underestimates the behavior 
in the tails as compared with measurement results. Since the 
calculation of detection probabilities for a communication sys¬ 
tem is primarily based on the tails of the pdf, underestimating 
this region significantly affects the accuracy of performance 
analysis. 

In an effort to address the shortcomings of the lognormal 
distribution, other statistical models have been further pro¬ 
posed to describe atmospheric turbulence channels under a 
wide range of turbulence conditions. These include the Nega¬ 
tive Exponential/Gamma model (also known widely as the K 
channel) (U, I-K distribution 10, log-normal Rician channel 
(also known as Beckman) ca, Gamma-Gamma ifTTl . M 
distribution ifT^ and Double-Weibull ifO . Particularly worth 
mentioning is the Gamma-Gamma model ifTll . ifT^ which has 
been widely used in the literature for the performance analysis 
of ESO systems, see e.g., ca, ca, along with the log-normal 
model. This model builds upon a two-parameter distribution 
and considers irradiance fiuctuations as the product of small- 
scale and large-scale fiuctuations, where both are governed by 
independent gamma distributions. In a more recent work by 
Chatzidiamantis et al in ina, the Double-Weibull distribution 
was proposed as a new model for atmospheric turbulence 
channels. Similar to the Gamma-Gamma model, it is based 
on the theory of doubly stochastic scintillation and considers 
irradiance fiuctuations as the product of small-scale and large- 
scale fiuctuations which are both Weibull distributed. It is 
shown in d that Double-Weibull is more accurate than the 
Gamma-Gamma particularly for the cases of moderate and 
strong turbulence. 

In this paper, we propose a new and unifying statistical 
model, named Double Generalized Gamma (Double GG), for 
the irradiance fiuctuations. The proposed model is valid under 
all range of turbulence conditions (weak to strong) and con- 
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tains most of the existing statistical models for the irradiance 
fluctuations in the literature as special cases. Furthermore, 
we provide comparison of the proposed model with Gamma- 
Gamma and Double-Weibull models. For this purpose, we 
use the set of simulation data from ifTTl . lIT^ for plane and 
spherical wave£]. Our model demonstrates an excellent match 
to the simulation data and is clearly superior over the other 
two models which show discrepancy from the simulation data 
in some cases. In the second part of the paper, we use this 
new channel model to derive closed form expressions for the 
BER and the outage probability of single-input single-output 
(SISO) and single-input multiple-output (SIMO) FSO systems 
with intensity modulation and direct detection (IM/DD). Our 
performance results can be seen as a generalization of the 
results in ll20ll - ll^ . 

The rest of the paper is organized as follows: In Section II, 
we propose Double GG distribution to characterize turbulence- 
induced fading. In Section III, we confirm the accuracy of our 
model through comparisons with simulation data for plane and 
spherical waves. In Section IV, we present the derivation of bit 
error rate (BER) and outage probability for SISO FSO system 
over Double GG channel. In Section V, we present BER 
expressions for ESO links with multiple receiver apertures. 
Einally, Section VI concludes the paper. 


II. Double GG Distribution 


The irradiance of the received optical wave can be modeled 
as im, ca / = Uy, where Ix and ly are statistically 
independent random processes arising respectively from large- 
scale and small scale turbulent eddies. We assume that both 
large-scale and small-scale irradiance fiuctuations are gov¬ 
erned by Generalized Gamma (GG) distributions ||24l Eq. (1)]. 
The pdfs of Ix ^ GG ( 71 , mi, f^i) and ly ^ GG ( 72 , m 2 , ^^ 2 ) 
are given by 


//. (4) = 

fly (ly) = 
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where 7 i > 0 , m^ > 0.5 and i = 1,2 are the GG 
parameters. The pdf of I can be written as 


fi(I) = I fiAI\Iy)fiy(Iy)dIy (3) 
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where //^ {I\Iy) is obtained as 


fi. 414) 


Iy(ni/mi 




^The simulation data in El, El was obtained through phase screen ap¬ 
proach which consists of approximating a three-dimensional random medium 
as a collection of equally spaced, two-dimensional, random phase screens that 
are transverse to the direction of wave propagation. In El, it was discussed 
in detail that such a numerical simulation approach contains all the essential 
physics for accurately predicting the pdf of irradiance (or, equivalently, log¬ 
normal irradiance), and shown that the simulation results provide an excellent 
match to known experimental measurements reported in El for both plane 
and spherical waves. The same set of simulation data was also used in El 
and El which respectively introduced Gamma-Gamma and Double-Weibull 
distributions as turbulence channel models. 


The integration in ([3]) yields 
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where [.] is the Meijers G-functiord defined in (251 
Eq.(9.301)], p and q are positive integer numbers that satisfy 
p/q = 71/72 and A{j;x) = x/j ,..., (x + j - l)/j . 
We name this new distribution as Double GG. Employing 
(261 Eq. (10)] and after some simplifications, the cumulative 
distribution function (cdf) of Double GG distribution can be 
obtained as 
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The distribution parameters 7 ^ and Qi i = 1,2 of the Double 
GG model can be identified using the first and second order 
moments of small and large scale irradiance fiuctuations. The 
latter are directly tied to the atmospheric parameters. Without 
loss of generality, we assume E (I) = 1 or equivalently 
E (Ix) = 1 and E (ly) = 1 . The second moment of irradiance 
is expressed as 


E (4) = E (4) E (4) = (1 + 4) (1 + 4) (7) 

where and are respectively normalized variances of Ix 
and ly. The n* moment of Ix (similarly/^) is given by 


E(i:) 


r (mi + n/ 71 ) 
mij r(mi) 


( 8 ) 


Therefore, by inserting the second order moment obtained 
from dSl) in (71), and considering that E (I) = 1, we have 


crl = 


^1 = 


a: = 


r (mi + 2/71 ) r (mi) 
4 (mi + 1/71 ) 
r (m 2 + 2/72 ) r (m 2 ) 
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where rrii is a distribution shaping parameter and found using 
curve fitting on the simulated/measured channel data. Note 
that in (951) and (EB, the variances of small- and large-scale 
uctuations (i.e., and cTy) are directly tied to the atmospheric 
conditions. Particularly, assuming a plane wave when inner 
scale effects are considered, the variances for the large- and 
the small-scale scintillations are given by (H Eqs. 9.46 and 


^Meijers G-function is a standard built-in function in mathematical 
software packages such as MATLAB, MAPLE and MATHEMATIC A. If 
required, this function can be also expressed in terms of the generalized 
hypergeometric functions using EH Eqs.(9.303-304)]. 
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(lnl-<lnl>)/a 


Fig. 1. Pdfs of the scaled log-irradiance for a plane wave assuming weak 
irradiance fluctuations. 



(Inl-<lnl>)/cj 


Fig. 2. Pdfs of the scaled log-irradiance for a plane wave assuming moderate 
irradiance fluctuations. 
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In ©, (Iq/Rq) is defined as 
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In the presence of a nite inner scale, the small-scale scintil¬ 
lation is again described by SlAi and the large-scale variance 
is given by ||4] Eq. 78] 


cr^ = exp 


0.04/321 


8.567/i 


7/6 


,8.56 + 7/, + 0.195/327/7/ 


where rji = 10.S9 (Rq/Iq), and Rq/Iq denotes the ratio of 
Fresnel zone to finite inner scale. 

For spherical waves in the absence of inner scale, and 
ay are given by [HI Eqs. 9.63 and 9.70] 
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Therefore, the parameters of the Double GG distribution are 
readily deduced from these expressions using only values 
of the refractive index structure parameter and inner scale 
according to the atmospheric conditions. The scintillation 
index can be further calculated as 


2 E (/ 2 ) 

<^1 = - 1 = (1 


E{IY 


\) (1 




(18) 


where /^q is the Rytov scintillation index of a spherical wave 
given by 


(15) 


We should emphasize that this distribution is very generic 
since it includes some commonly used fading models as 
special cases. For 7 ^ ^ 0, ^ cx). Double GG pdf coincides 

with the log-normal pdf. For 7 ^ = 1, (2^ = 1, it reduces 
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(lnl-<lnl>)/CT 

Fig. 3. Pdfs of the scaled log-irradiance for a plane wave assuming strong 
irradiance fluctuations. 


TABLE I 

NRMSE FOR DIFFERENT STATISTICAL MODELS AND 
TURBULENCE CONDITIONS DEEINED IN FiGS. 1-6 



Gamma-Gamma 1111 

Double-Weibull 1131 

Double GG (Proposed) 

Fig. 1 

2% 

7% 

0.6% 

Fig. 2 

2.8% 

1.2% 

0.8% 

Fig. 3 

1.2% 

1% 

0.8% 

Fig. 4 

0.3% 

10% 

0.3% 

Fig. 5 

19% 

8.7% 

1.5% 

Fig. 6 

4.8% 

2.4% 

1.7% 


to Gamma-Gamma while for rrii = 1, it becomes Double- 
Weibull. For 7 ^ = 1, = 1, m 2 = 1, it coincides with the K 

channel. 

III. Verification of the Proposed Channel Model 

In this section, we compare the Double GG distribution 
model with simulation data of plane and spherical waves pro¬ 
vided respectively in ini and ca. In ini, Flatte et al carried 
out exhaustive numerical simulations and published the results 
assuming plane wave propagation through homogeneous and 
isotropic Kolmogorov turbulence. In ifTSl . Hill and Frehlich 
presented the simulation data for the propagation of a spher¬ 
ical wave through homogeneous and isotropic atmospheric 
turbulence. The turbulence severity is characterized by Rytov 
variance (cr^yioy) which is proportional to the scintillation index 
IQ- We emphasize that the same data set was also employed 
in na and Ha which have introduced the Gamma-Gamma 
and Double-Weibull fading models. 

A. Plane Wave 

Figs. 1-3 compare the Gamma-Gamma, Double-Weibull and 
Double GG models under a wide range of turbulence condi¬ 
tions (weak to strong) assuming plane wave propagation. In 
these figures, the vertical axis of the figure represents the log- 
irradiance pdf multiplied by the square root of variance. The 
logarithm of irradiance was particularly chosen to illustrate 



(lnl-0.5a)/a 

Fig. 4. Pdfs of the scaled log-irradiance for a spherical wave assuming weak 
irradiance fluctuations. 

the high and low irradiance tails ifTTl . Thus, sensitivity to the 
small irradiance fades is increased, while sensitivity to large 
irradiance peaks is decreased. The pdf plots were also scaled 
by subtracting the mean value to center all distributions on 
zero and dividing by the square root of variance. 

In Fig. 1, we assume weak turbulence conditions which are 
characterized by = 0.1 and Iq/Rq = 0.5. The values 

of the variances of small and large scale fluctuations, (cr^ and 
a^) are calculated from (HHi and (Ha. Using (iiii, (ibii and 
(flOl) . the Double-GG parameters are obtained as 71 = 2.1, 
72 = 2.1, mi = 4, m 2 = 4.5, Oi = 1.0676 and Q 2 = 1-06 
where p = q = 1. We further employ normalized root-mean- 
square error (NRMSE) as a statistical goodness of fit test. 
Table I provides the NRMSE results for different statistical 
models. According to Table I and Fig. 1, both Gamma-Gamma 
and Double-Weibull distributions fail to match the simulation 
data. On the other hand, the proposed Double GG distribution 
follows closely the simulation data. 

In Fig.2, we assume moderate irradiance fluctuations which 
are characterized by cr^y^Q^ = 2 and Iq/Rq = 0.5. The 
parameters of the Double GG distribution for this case are 
obtained as 71 = 2.1690, 72 = 0.8530, mi = 0.55, m 2 = 
2.35, = 1.5793 and U 2 = 0.9671. In the calculations, 

p and q are chosen as p = 28 and q = 11 in order to 
satisfy p/q = 71 / 72 . Among the three distributions under 
consideration, the proposed Double GG model provides the 
best fit to the simulation data. It is apparent that Gamma- 
Gamma fails to match the simulation data particularly in the 
tails. As Table I demonstrates, the accuracy of Double Weibull 
is better than that of Gamma-Gamma, but slightly inferior to 
our proposed distribution. 

In Fig. 3, we assume strong irradiance fluctuations which are 
characterized by cr^y^^yQ = 25 and Iq/Ro = 1- The parameters 
of the Double GG distribution are calculated as 71 = 1.8621, 
72 = 0.7638, mi = 0.5, m 2 = 1.8, Qi = 1.5074 and U 2 = 
0.9280 where p and q are chosen as 17 and 7 respectively. 































5 



Fig. 5. Pdfs of the scaled log-irradiance for a spherical wave assuming 
moderate irradiance fluctuations. 



Fig. 6. Pdfs of the scaled log-irradiance for a spherical wave assuming strong 
irradiance fluctuations. 

The Double GG model again provides an excellent match to 
the simulation data and as it is clear from Table I, its accuracy 
is better than Gamma-Gamma and Double Weibull. 

B. Spherical Wave 

Figs 4-6 compare the Gamma-Gamma, Double-Weibull 
and Double GG models under weak, moderate and strong 
turbulence conditions assuming spherical wave propagation. 
These pdfs are plotted as a function of (In / + 0.5cr^)/cr 
where a is the square root of the variance of In /. The y-axis 
depicts the log-irradiance pdf multiplied by a. 

In Fig. 4, we consider spherical wave propagation and 
assume weak turbulence which are characterized by = 

0.06 and Iq/Rq = 0. The parameters of Double GG are evalu¬ 
ated using the variances of small and large scale fluctuations. 





(C) 



Fig. 7. Outage probability as a function of for a) Plane wave - 

^Rytov = 2, lo/Ro = 0.5, b) Plane wave-(T|„t„„ = 25, Iq/Rq = 1, c) 
Spherical wave - = 2, Iq/Rq = 0, d) Spherical wave - = 5, 

Iq/Rq = 1 





















































6 


(ay and cr^) for spherical waves. The values of these variances 
are given by (US , (O and dnii. Therefore, employing (9) 
and ([TOb we obtain mi = 34.24, m 2 = 32.79, 71 = 72 = 
= (^2 = 1 where p and q are equal to 1. It can be noted 
that in this case, the Double GG coincides with the Gamma- 
Gamma distribution. It is apparent that both Gamma-Gamma 
and Double GG distributions provide an excellent match to the 
simulation data while the Double-Weibull distribution fails to 
match the simulation data. 

In Fig. 5, we assume moderate irradiance fluctuations, which 
are characterized by = 2 and Iq/Rq = 0. The param¬ 

eters of the Double GG model for this case are calculated 
as 71 = 0.9135, 72 = 1.4385, mi = 2.65, m 2 = 0.85, 
Qi = 0.9836 and D 2 = 1.1745 where p and q are selected as 
7 and 11 respectively. It is clearly observed that the Double 
GG model provides a better fit with simulation data, especially 
for small irradiance values. 

In Fig. 6 , we assume strong irradiance fiuctuations which are 
characterized by = 5 and Iq/Rq = 1. The parameters 

of the Double GG model are calculated as 71 = 0.4205, 
72 = 0.6643, mi = 3.2, m 2 = 2.8, fli = 0.8336 and 
D 2 = 0.9224 where p and q are chosen as 7 and 11 
respectively. It is apparent from this figure and Table I that 
both Gamma-Gamma and Double-Weibull distributions fail to 
match the simulation data. On the other hand, the proposed 
Double GG distribution follows closely the simulation data. 

IV. Performance Evaluation 
A. Outage Probability Analysis of SISO FSO System 

Denote Rt as a targeted transmission rate and assume 
"fth = {Rt) as the corresponding signal-to-noise ratio 
(SNR) threshold in terms of the instantaneous channel capacity 
for a particular channel realization. Therefore, the outage 
probability is calculated by Pout {Rt) = Pr {y < yth) 1 ^ - 
|[28]| . If SNR exceeds yth, no outage happens and the receiver 
can decode the signal with arbitrarily low error probability. For 
the system under consideration, the instantaneous electrical 
SNR can be defined as 7 = {r]I)‘^/N q, while the average 
electrical SNR is obtained as 7 = rf/Nn since E {!) = 1. 
Therefore, / can be expressed as / = sjyjy. After the 
transformation of the random variable, /, the cdf of 7 can 
be easily derived from ^ and setting 7 = jth therein, we 
obtain the outage probability as in ([T9l) at the top of the next 
page. 

(US can be seen as a generalization of earlier outage 
analysis results in the literature. Specifically, if we insert 
7 i = 1 and = 1 in ([19]) , we obtain the outage probability 
expression reported in ||20| Eq* (15)] for Gamma-Gamma 
channel. Setting = 1 in ([19]), we recover Eq (16) of |[T3]| 
derived for Double-Weibull channel. Similarly, for 7 ^ = 1 , 
Oi = 1 and m 2 = 1, (O reduces to (3) of reported for 
the K channel. 

Based on the derived expression in ([T9]) . Figs. 7.a-d present 
the outage probabilities of a SISO FSO system for different 
degrees of turbulence severity. We adopt the same parameters 
used in Figs 2-3 and 5-6 and consider the following four 
cases: a) Plane wave with cr^y^oy = 2 and Iq/Rq = 0.5, b) 





(c) 



Fig. 8 . Average BER as a function of 7 a) Plane wave - cr‘^ytov — 
Iq/Rq = 0.5, b) Plane wave-a^^^^^ = 25, Iq/Rq = 1, c) Spherical wave - 
(T%ytov — 2, Iq/Rq = 0 , d) Spherical wave - cr‘^ytov — ^o/Ro = 1 
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{n2pf{^iqy ^ (Q : mi), a (p : m 2 ), 0 


(19) 


Plane wave with cr^ytoy = 25 and Iq/Rq = 1, c) Spherical 
wave with cr^y^Qy = 2 and Iq/Rq = 0, d) Spherical wave with 
^Rytov = ^ /o/^o = 1 . It is observed from Fig. 7.a 

that an SNR of 37.8 dB is required to achieve a targeted 
outage probability of 10“^. As the turbulence strength gets 
stronger (see Fig. 7.b), the required SNR to maintain the same 
performance climbs up to 50.5 dB. Similarly, for spherical 
waves, SNRs of 36.8 dB and 50.9 dB are respectively required 
for moderate and strong turbulence conditions. In these figures, 
we further include the outage results for Double-Weibull and 
Gamma-Gamma for comparison purposes. As expected from 
the earlier comparisons of their pdfs, the outage performance 
over Double-Weibull and Double-GG for plane wave (See 
Figs. 7.a and 7.b) are similar while the Gamma-Gamma model 
overestimates the outage performance. On the other hand, 
the superiority of Double GG is more obvious for spherical 
wave (See Figs 7.c and 7.d), particularly for strong turbulence 
conditions, where the outage performance plots of Double- 
Weibull and Gamma Gamma significantly deviate. 


B. BER Analysis of SISO FSO System 

In this section, we present the BER performance analysis of 
an FSO system with on-off keying (OOK) over the proposed 
Double GG channel. The received optical signal is written as 

y = ylx -h n (20) 


where x represents the information bits and can be either 0 
or 1, n is the Additive White Gaussian noise (AWGN) term 
with zero mean and variance = Aq/2 , 7 / is the optical-to- 
electrical conversion coefficient and / is the normalized irra- 
diance whose pdf follows ©. Conditioned on the irradiance, 
the instantaneous BER for OOK is given by ll22l 

P„..=0^5erfc(^) , 21 ) 

where erfc (.) stands for the complementary error function 
defined as 

erfc(x) = —= / e~^ dt. ( 22 ) 

V ^ J X 

The average BER can be then calculated by averaging (121]) 
over the distribution of /, i.e.. 


Pe = 



0.5 erfc 



dl 


(23) 


The above integral can be evaluated in closed form by ex¬ 
pressing the erfc (.) integrand via a Meijers G-function using 
(301 Eq. (8.4.14.2)], O Eq. (8.2.2.14)] and (lH Eq. (21)]. 
Thus, a closed-form solution is obtained as in (l24l) at the top 
of the next page. 


In ([24]) k and I are positive integer numbers that satisfy 
P 72/2 = l/k and {y,x) is defined as 


iy,x) 

= A A 


(25) 


y-x 
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.Ak, 


y-l-x 


,aU, 


1 — X 


The derived BER expression in (|24]) can be seen as a 
generalization of earlier BER results in the literature. If we 
insert 7 ^ = 1 and = 1 in (|24]) . we obtain the BER 
expression derived in (211 Eq. (9)] under the assumption of 
Gamma-Gamma channel. Setting = 1 in (1241) . we obtain 
Eq (15) of f\3\ derived for Double-Weibull channel. On the 
other hand, for 7 ^ = 1, = 1 and m 2 = 1, (1241) reduces to 

(12) of reported for the K-channel. 

In an effort to have some further insights into system per¬ 
formance, we investigate the asymptotical BER performance 
in the following. Eor large SNR values, the asymptotic BER 
behavior is dominated by the behavior of the pdf near the 
origin, i.e. // (/) at / ^ 0 [321. Thus, employing series 
expansion corresponding to the Meijers G-function (331 Eq. 
(07.34.06.0006.01)], the Double GG distribution given in (0) 
can be approximated by a single polynomial term as 
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In ([26]) . bk and hj are defined as 
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Therefore, based on (| 2 l, the average BER can be well 
approximated by 
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Erom (|30l(, it can be readily deduced that the diversity order 
of SISO ESO system is given by 0 . 5 p 72 mm{mi/q ^m^/p). 
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Fig. 9. Comparison of the average BER between SISO and SIMO with 
optimal combing for plane wave as defined in Fig. 3. 



Fig. 10. Comparison of the average BER between SISO and SIMO with 
optimal combing for spherical wave as defined in Fig. 5. 


It is observed from Fig.8.a that an SNR of 51.1 dB is 
required to achieve a BER of 10“^ for a plane wave in moder¬ 
ate turbulence conditions. For stronger turbulence conditions, 
the required SNR to achieve the same BER performance is 
68.2 dB as seen from Fig. 8.b. For spherical waves, SNRs of 
49.8 dB and 63.8 are respectively required for moderate and 
strong turbulence conditions. Comparison with the expressions 
presented for other channel models reveals that the Gamma- 
Gamma model significantly overestimates the performance. 
Similar to earlier observations on the outage analysis, the 
superiority of Double GG is more obvious for spherical wave. 
As observed from Figs 8.c and 8.d, the performance plots 
of Double-Weibull and Gamma Gamma considerably diverge 
particularly for strong turbulence conditions. Furthermore, it 
can be clearly seen that the asymptotic results are in excellent 
agreement with exact analytical results within a wide range 
of SNR showing the accuracy and usefulness of the derived 
asymptotic expression given in (1^ . 


V. FSO Links with Receive Diversity 


The optimum decision metric for OOK is given by 1^ 

on 

P(r|on,I„)^P(r|off,I„) (31) 

off 


where r = (ri,r 2 ,tat) is the received signal vector and 
is the fading channel coefcient which models the channel from 
the transmit aperture to the n* receive aperture. Following 
the same approach as Ea, ED the conditional bit error 
probabilities are given by 


Pe(0ff |4)=Pe(0n 14) 



(32) 


The average BER can be then obtained as 


4imo,oc ^ J ^ 



(33) 


As it can be noticed from Section IV, the performance of 
a SISO FSO link over moderate and strong atmospheric tur¬ 
bulence is quite poor. To address this issue, multiple transmit 
and/or receive apertures can be employed and the performance 
can be improved via diversity gains. In the following, we 
assume that multiple receive apertures are employed and 
present the BER derivations under the assumption that optimal 
gain combining is used. 


where fi (I) is the joint pdf of vector I = (p, / 2 ,..., /at). 
The factor N in ([3^ is used to ensure that sum of the N 
receive aperture areas is the same as the area of the receive 
aperture of the SISO link. 

Eq. ([3^ does not yield a closed-form solution and requires 
N-dimensional integration. Nevertheless, the Q-function can 
be well-approximated as Q{x) ^ e“^/12 -f- e“^/4 ll34ll 
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(34) 


and thus the average BER can be obtained as 
-PsiMO.OC ~ n / ( 4/V'^” ) 

+ i n /°° (^^n) din (35) 

The above integral can be evaluated by first expressing the ex¬ 
ponential function in terms of the Meijer G-function presented 
in EH eq. (11)] as 

exp (-a;) = GJ;? [a; |o ] (36) 

Then, a closed-form expression for (l35ll is obtained using ll3Tl 
Eq. (21)] as 


1 - - 1 , , 

-PsiMO,oc ~ f[ A (n, 4) -I- - P[ A (n, 3) (37) 

n=l n=l 

where A (n, v) is defined in (l34li at the top of the page. 

In El, In and kn are positive integer numbers that satisfy 
Pn72,n/2 = In/kn , and an and cUn n e {1,2,..., N} are 
defined as 


an 

UJn 


m 2 

j 2 ,nPn 


t, + 1/2 


mi^n 1/2 
Qn 


(277) 


l-(Pn+<7rr)/2 


r ( 7711 ^ 72 ) r (Tfl2,n) 


(38) 

(39) 


The derived expression in (ITtI) includes the previously reported 
result in Ea for K channel as a special case. 

Based on the approximation in (l26l) . the corresponding 
closed-form asymptotic solution for (1351) can be obtained as 


\ 1 

^SIMO,OC_asy ~ — n -^asy («,4)+- n -^asy {n,3) (40) 

n=l n=l 

where Aasy (ri, 11 ) is defined as 


Pn+Qn 

-^asy 11) — T (^jf,n ^k,n) 

i=i 


^ r (Pn 72 ,n^/c,n) 


Pn 72 ,nbk,n 





(41) 


Therefore, the diversity order of FSO links with A/ receive 
apertures employing optimal gain combining is obtained as 

N 

0.5 y] Pn 72 ,nniin {mp„/g„,m 2 ,„/p„}. 

n=l 

Figs. 9-10 illustrate the BER performance of the SIMO 
FSO system under consideration. We present approximate 
analytical results which have been obtained through ([JTl) and 
go]) along with the Monte-Carlo simulation of (1^ . As clearly 
seen from Figs. 9-10, our approximate expressions provide an 
excellent match to simulation results. As a benchmark, the 


BER of SISO FSO link is also included in these figures. 
It is observed that receive diversity signicantly improve the 
performance. For instance, at a target bit error rate of 10“^, 
we observe performance improvements of 26.8 dB and 39.6 
dB respectively for with N = 2 and 3 with respect to the SISO 
transmission over Double GG turbulence channels defined in 
Fig.3. Similarly, for Double GG channels defined in Fig. 5, 
at a BER of 10“^, performance improvements of 19 dB and 
25.1 dB are achieved for SIMO links with N = 2 and 3 
compared to the SISO deployment. It can be further observed 
that asymptotic bounds on the BER become tighter at high 
enough SNRs confirming the accuracy and usefulness of the 
asymptotic expression given in (l40b . 

VI. Conclusions 

In this paper, we have introduced a new channel model, 
so called Double GG, which accurately describes irradiance 
fiuctuations over atmospheric channels under a wide range 
of turbulence conditions. It is based on the theory of doubly 
stochastic scintillation and considers irradiance fiuctuations as 
the product of small-scale and large-scale fiuctuations which 
are both Generalized Gamma distributed. We have obtained 
closed-form expressions for the pdf and cdf in terms of 
Meijers G-function. Comparisons with the Gamma Gamma 
and Double-Weibull have shown that the new model pro¬ 
vides an accurate fit with numerical simulation data for both 
plane and spherical waves. Using the new channel model, 
we have obtained closed-form expressions for the BER and 
the outage probability of SISO and SIMO FSO systems. We 
have demonstrated that our derived expressions cover many 
existing results in the literature earlier reported for Gamma- 
Gamma, Double-Weibull and K channels as special cases. 
Based on the asymptotical performance analysis, we have 
further derived diversity gains for SISO and SIMO FSO 
systems under consideration. 
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